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Abstract 

The structure of polymer coils near interfaces between coexisting phases of sym- 
metrical polymer mixtures (AB) is discussed, as well as the structure of symmetric 
diblock copolymers of the same chain length iV adsorbed at the interface. The prob- 
lem is studied by Monte Carlo simulations of the bond fluctuation model on the 
simple cubic lattice, choosing N = 32 and lattice linear dimensions L x D x L up 
to 512 x 64 x 512, and using massively parallel computers (CRAY T3D). While ho- 
mopolymer coils in the strong segregation limit are oriented parallel to the interface, 
the diblocks form "dumbbells" oriented perpendicular to the interface. However, in 
the dilute case ("mushroom regime" rather than "brush regime"), the diblocks are 
only weakly stretched. Distribution functions for monomers at the chain ends and 
in the center of the polymer are obtained, and a comparison to the self consistent 
field theory is made. 



1 Introduction 



Blending of polymers is useful to obtain materials with improved properties [1]. 
However, most pair (A,B) of chemically different homopolymers do not mix: 
any unfavorable enthalpy per monomer, multiplied by the large number of 
monomers in a chain ("chain length" N^, Nb] for simplicity only " symmetric 
mixtures", Na = Nb = N are treated here), will exceed the entropy of mixing 
(which is of order fc^T, T being the absolute temperature) [2] . Therefore the 
widely used polymer blends are not homogeneous on mesoscopic scales, rather 
they are fine dispersion of one polymer in another [3], i.e. the material is full 
of AB-interfaces, and the detailed structure of the polymers in the interfacial 
region clearly have a crucial importance for the understanding the application 
properties of the blend. 

In practice one often uses block copolymers [4,5] as polymeric surfactants that 
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control the size of the minority droplets in the dispersion. The classical expla- 
nation[l,6] for the usefulness of the block copolymer to produce fine enough 
dispersion is based on the fact that the adsorption of block copolymers at 
the AB interface lowers the interfacial tension[7,8]. However, this mechanism 
was recently questioned[3,9,10], and it was argued the primary effect of the 
adsorbed block copolymer is a kinetic one, i.e. they prevent the small droplets 
to a large extent from coagulation. 

For a more complete understanding one clearly needs detailed knowledge of in- 
terfacial properties, both with and without adsorbed copolymers, including the 
detailed configurations of the chains when they are confined to the interfacial 
region. While numerous elegant experiments have yielded already interesting 
information on this problem (e.g. [11-16]), computer simulation [17-1 9] is a 
tool that can provide a much more detailed and complete picture. Although 
the models accessible to simulation [20] are very idealized, they are well defined 
and all the quantities needed to stringently test corresponding theories such 
as the self consistent field (SCF) theory[8,18,19,21-23] can easily be estimated 
from such simulations. There is no problem with fitting of phenomenological 
parameters to experiments. In the present paper, we focus on the configura- 
tional properties of the chains in the interfacial region of a strongly segregated 
mixture, considering also adsorbed diblock copolymers of the same chemical 
nature AB, symmetric composition (/ = 1/2) [4,5], and the same chain length 
as the homopolymers. 



2 Some comments on the model and the simulation method 

While simulational efficiency requires to sacrifice chemical detail [20], there is 
ample evidence [5, 20] that a lot of useful insight can be gained already from 
coarse-grained models of flexible polymers such as the bond-fluctuation model [24- 
26]. Each effective monomer blocks a cube of 8 neighboring sites from further 
occupancy on a simple cubic lattice. Effective monomers are connected by ef- 
fective bond vectors of length 2, \/5, VQ, 3 or y/lO in units of the lattice spacing 
a = 1. (Each effective bond represents a group of n « 3 — 5 subsequent C — C- 
bonds along the backbone of the chain.) Using a chain length N = 32, as done 
here, corresponds to a degree of polymerization of 100 — 160 in a real poly- 
mer. For a volume fraction $ = 0.5 of occupied sites, properties of a dense 
melt are well reproduced [25, 26]. Interactions between monomers are modeled 
by energy parameters €aa = £bb = —£ab = —ksT, if the distance between 
monomers does not exceed \/6- The statistical segment length b in the rela- 
tion for the radius of gyration R g = b^jN/6 is b = 3.05 (R g ~ 7 for N = 32). 
The Flory-Huggins parameter is x = ^ z eff^ where z e ff ~ 2.65 denotes the 
effective coordination number in the bulk[17]. The normalized compressibility 
is k B Tna^ = 4.1 [19]. The bulk phase diagram[25], the correlation length of 
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Fig. 1. Monomer density normalized by the bulk density p m (y)/ ' Pbulk plotted vs. 
y/wssL where the (theoretical [21]) interfacial width in the strong segregation limit 
(SSL) is wssl = b/ (6X) 1 / 2 ~ 1.71 at e = 0.1. Shown are the total density p, A and B 
monomers separately (pa,Pb), homopolymers irrespective of their nature (pn), and 
A and B monomers belonging to a copolymer (pAd Pbc)- The monomer profiles for 
a pure homopolymer system without copolymers are p A ,p° B . From Ref.[19]. 

concentration fluctuations[27], etc. have been determined in previous work. 
Interfaces are studied in a L x D x L geometry where lattice linear dimensions 
L parallel to the interface are L = 512, and perpendicular D = 64. Choosing 
boundary conditions periodic in x,z directions and "antiperiodic" in y direction 
(i.e. an A-chain part leaving the box at y = —D/2 reenters at y = +D/2 + 1 
as a B-chain, etc. ) one maintains a system with a single interface at y — 0. 
The system is initialized (for homopolymer interfaces) by choosing chains that 
have their center of gravity in the left part of the box (y < 0) as B-chains, 
and in the right part as A-chains. This corresponds to complete segregation 
between A and B in the bulk, appropriate for large enough e (e.g. e = 0.1). 
For weakly segregated mixtures, rapid equilibration of bulk concentrations 
PaiPb is achieved by use of semi-grandcanonical algorithms [17, 25]. However, 
very long runs to equilibrate the model system are nevertheless required, in 
order to allow for a build-up of long range capillary wave fluctuations of the 
interface. Note, that the interfacial position itself can fluctuate and therefore 
the origin of the y-axis is always fixed at the center of the (instantaneous) 
interfacial profile for all "measurements". 

For simulating interfaces with adsorbed copolymers, we initialize the system 
using an equilibrated homopolymer interface. Then, by choosing 1024 chains 
(out of 32768) whose center of mass is in the interval [—5, +8] with 8 = 3 or 
9, we "transform" them into copolymers. The configurations are equilibrated 
with 2.5 ■ 10 5 attempted moves per monomer, using a random hopping algo- 
rithm^]. For e = 0.1 (where the concentration of the copolymer in the bulk 
is 0.04% [28]) no effect of varying 8 was found. We average over 86 configura- 
tions, taken every 10 4 attempted moves per monomer. We estimate that under 
the chosen conditions about 30% of the interfacial area is covered by copoly- 
mers, i.e. we study the dilute case ("mushroom regime" rather the "brush 
regime" [7]). 
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Fig. 2. Homopolymer segment density profiles p s (y) plotted vs. y/wssL, for 
monomers in the middle of the chain (pi/2), at the chain ends (p e ), and at one 
and three quarters of the chains (P1/4) and of all homopolymer monomers (pn)- 
Also shown is the distribution of midpoints between the two ends of the homopoly- 
mer p ee . Units are the bulk density pj, or p s ^ = , as indicated. From Ref.[19] 



Fig. 3. Block copolymer segment density profiles for densities of A,B monomers in 
the middle of the chains (pi/2Aj P\I2b)i at the chain ends (p e A, PeB), at one and three 
quarters of the chain (pi/4,4, Pi/ab)-, and of all copolymer monomers (pcA, Pcb)- The 
inset shows the results of the SCF theory for A monomers. From Ref. [19]. 

3 Density profiles of monomers and single segments 



Fig. 1 shows profiles of A and B monomer densities in systems with and with- 
out copolymers. It is seen that in this dilute case A and B monomer profiles 
are the same, whether copolymers are present or not. The slight dip of the 
total density, i.e. the enrichment of "free volume" in the center of the inter- 
face, can be explained quantitatively in terms of the nonzero compressibility 
of the model[17]. Qualitatively, the distribution of the copolymer monomers 
agree with experimental results, monomers of type A being concentrated in 
the A-rich phase, and monomers of type B in the B-rich phase. 
A more detailed picture results from the distribution of single chain segments 
(Figs. 2,3). Chain ends of homopolymers are enriched at the interface rela- 
tive to the total homopolymer density profile, as in the case without copoly- 
mers[17]. Profiles of inner segments (pi/2,Pi/4) are close to the total density 
profile. However the "geometrical midpoint" distribution p ee has a minimum 
at the interface and a maximum a gyration radius away: this tells that typi- 
cally chains have only one end at the interface and not both ends! 
As expected, copolymer middle segments concentrate at the interface, whereas 
chain ends stretch out in their favorite bulk phase. The distribution of monomers 
in the middle of a block resembles that of the total monomer distribution of 
the corresponding species of the copolymer. All these features are in qualita- 
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Fig. 4. Orientational asymmetry parameter q e (y) of the end-to-end distance plotted 
vs. y at 6 values of e as indicated. From Ref.[17]. 

tive agreement with SCF theory (the underestimation of the interfacial width 
by the SCF theory causes some minor qualitative deviations, of course). 



4 Orientation of chains and of bonds 

Although homopolymers typically have only one chain end right near the inter- 
facial center, there is nevertheless a pronounced effect on the orientation of the 
chains. This is seen by studying the orientational order parameter (Fig. A)q e 
for the end-to-end distance, q e (y) = [3(Ry) y — {R 2 )y}/2(R 2 ) y where the outer 
index y at the brackets (• • -) y denotes the center of gravity y-coordinate of the 
considered chain, the inner indices x, y, z denote the Cartesian components 
R = (R x ,Ry,R z ). One sees that q e (y) ~ (i.e. random orientation) in the 
weak segregation limit (note that the critical point occurs for e c = 0.0143(25]), 
while the strongly negative values of q e (y ~ 0) in the SSL indicate that paral- 
lel components of R are much larger than perpendicular ones. A similar effect 
occurs for the gyration radius as well, but there is hardly any effect for in- 
dividual bond vectors[17]. Analyzing the mean square gyration components 
in the eigensystem of the gyration tensor, one easily finds that Fig. 4 is not 
primarily an effect of coil deformation (there occurs only a slight shrinkage of 
chain linear dimensions) but due to coil orientation (the instantaneous soap- 
shaped coils are oriented with their two longer axis more or less parallel to 
the interface). 

Copolymers show the inverse behavior, they stretch in the direction per- 
pendicular to the interface (Fig. 5). Both SCF and Monte Carlo calculations 
predict that the effect is strongest for the copolymers centered at about one 
to two radii of gyration away from the interface, and much weaker for those 
chains located in the wings of the concentrations profile or at the middle of 
the interface. Thus one can picture the copolymer in the latter case as con- 
sisting of two almost independent blocks, which hardly feel the effect of being 
linked together. One finds that the vectors connecting the end of the single 
A or B blocks are on average hardly oriented[19]. Only the blocks centered 
deep in their majority phase (y/wssL ~ ±7) stretch perpendicular to the 
interface, since they are pulled toward the interface by the other copolymer 
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Fig. 5. Mean square end-to-end vector components (R 2 ) with i = x, y, z in units of 
the average bulk value b 2 N/3, plotted vs. the distance of the center of the end-to-end 
vector from the interface y, at e = 0.1. Results are shown for homopolymers and 
copolymers and compared to the SCF prediction for the latter. From Ref.[19]. 



Fig. 6. a) Orientational order parameter q for various bonds of copolymers as 
function of their position y. Included are end bonds, link bonds (from 16th to 
17th monomer), bonds next to the link bonds, and bonds in the middle of a block 
(from 8th to 9th and from 24th to 25th monomer) b) Same as a) but from a SCF 
calculation. From Ref.[19]. 

end. Also the vector D connecting the centers of the mass of the A and B 
blocks is on average oriented and strongly stretched in the negative y direc- 
tion((4) = -3.8,v^>=4.5). 

Thus in this dilute situation single blocks are mostly not oriented at all: the 
perpendicular orientation of whole copolymers results from the arrangement 
of the two constituent blocks. Diblock copolymers resemble dumbbells consist- 
ing of two mildly perturbed homopolymer coils, and thus have a conformation 
similar to that found in the copolymer melts in the disordered phase near the 
order-disorder transition [29]. 

This picture is corroborated by the orientation of individual bonds, q(y) = 
[3(by)y — {b 2 )y}/2{b 2 ) y , see Fig. 6. For homopolymers q(y) stays extremely small 
throughout, q(y ~ 0) ~ —0.01 to —0.02 due to the on average parallel ori- 
entation of coils close to the interface. One finds a similar behavior for the 
end bonds or the bonds in the block middle for the copolymers, while the link 
bonds are oriented perpendicular to the interface (q(y) > 0, in particular if y 
is one or two gyration radii away from the interface center). This behavior is 
reproduced by the SCF theory only qualitatively. Note that the SCF profiles 
for q(y) clearly reflect the two different length scales: the width Wssl controls 
the extent of the central dip, while the gyration radius controls the overall 
width of the region with nonzero q. Due to capillary waves fluctuations, this 
distinction is smeared out in the Monte Carlo results. 
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5 Discussion 



As far as corresponding experimental results are available, e.g. from the work of 
Russell and coworkers [14], experiment and simulation agree qualitatively. On 
the other hand, by obtaining simultaneous information on both density profiles 
of various monomers along the chains and on bond orientations we clearly can 
go beyond experiment, resulting in a very clear and detailed "picture" how 
homopolymer and block copolymer configurations look like. But clearly the 
present work, restricting attention to a single chain length (chosen also the 
same for both homopolymers and copolymers) is a first step only, and one 
must already expect a different picture when one studies more concentrated 
block copolymer layers at interfaces. Also various asymmetries (different chain 
length of the blocks, different chain stiffnesses, etc. ) deserve attention. 
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Figure 5 
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